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Abstract. We study the properties of the Polyakov loop in strongly-coupled gauge 
plasmas that are conjectured to be dual to five dimensional theories of gravity coupled 
to a nontrivial single scalar field. We find a gravity dual that can describe the 
thermodynamic properties and also the expectation value of the Polyakov loop in 
the deconfincd phase of quenched SU(3) QCD up to 3T C . 



1. Introduction 

The expectation value of the Polyakov loop is an order parameter for the deconfinement 
phase transition in SU(N C ) gauge theories without fermions |1] since Z(N C ) symmetry 
is broken in the deconfined phase. At large N c the phase transition becomes a strong 
first-order phase transition [2] and the Polyakov loop jumps from zero to a finite value 
at T c (which is also known as the Gross- Witten point [3]). We define the path-ordered 
Polyakov loop as 

L(f) = J-T r Peifo^M^dr ^ 

where is the non-Abelian gauge field operator and the trace is over the fundamental 
representation of SU(N C ). According to the equation above, the loop can be seen as 
the trace of the propagator of an infinitely heavy quark in imaginary time. A less 
abstract way to understand this operator involves the connection between the loop and 
the free energy of a single, isolated (and infinitely massive) heavy quark in the medium 
® F Q (T), i.e., 

£=\(L)\ = e- F ^ T . (2) 

Thus, in the confined phase of SU(N C ) Yang-Mills (YM) theories the free energy diverges 
(because of Wilson's area law) and £ then vanishes. If dynamical fermions are added then 
I is small but nonzero in the hadronic phase because QQ pairs can be created thermally. 
The loop in QCD was computed in perturbative theory a long time ago and the leading 
correction to its asymptotic value (which is unity because of asymptotic freedom) is 
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positive and of order g 3 [5]. This implies that in perturbation theory I approaches unity 
from above. We will show in this paper that gauge theories with gravity duals that 
can mimic the current lattice predictions for the thermodynamics of pure glue can also 
describe £. 

Because of cluster decomposition, £ can also be obtained from the disconnected 
piece of the Polyakov loop correlator 



separated by a distance r in the medium (for a recent review of quarkonia at finite 
temperature see [H]). 

2. Gravity Duals and the Breaking of Conformal Invariance 

The gauge/string duality states that the complicated non-Abelian dynamics of four 
dimensional gauge theories can be understood in terms of the physics of higher 
dimensional theories of gravity [TJ |S]. This duality establishes that when the gauge 
theory is strongly-coupled the corresponding gravitational description is weakly-coupled, 
which makes this whole idea hard to falsify. However, there is in fact enough evidence 
that the main ideas behind this conjecture are right [9]. 

The most well studied example of this correspondence is the duality between 4- 
dimensional A* = 4 Supersymmetric Yang-Mills (SYM) and type IIB string theory on 
AdSs ® S5 [7\. In this case, the gauge theory at large N c is in the regime of strong 
t'Hooft coupling A ^> 1 and the typical "stringy" coupling in the gravitational theory 
is a /R 2 = 1/y/X < 1, where y/a' is the string length and R is the radius of AdS§. 
Here we refer the 5-dimensional spacetime in the gravity duals as the "bulk" and the 
4-dimensional flat spacetime where the gauge theory is defined as the "boundary". 
Note that the boundary of AdS$ is a flat 4d spacetime. Also, because in YM the only 
degrees of freedom are adjoint gluons, no other adjoint scalar fields in the gauge theory 
are included here and our working assumption is that the gravity dual theory may be 
properly defined in a 5 dimensions at least when it comes to the observables discussed 
in this paper. 

Gauge theories at finite temperature T in general correspond to dual geometries 
with a black brane and Hawking temperature T HW = T [TO]. A/" = 4 SYM is exactly 
conformal invariant in the vacuum and at finite temperature the ratio p/T 4 (where p is 
the pressure) does not depend on T. Here we consider the work of Gubser and Nellore 
in [H] that describes gravity duals where conformal invariance is broken in the infrared 
but restored when T — > 00. This is described using a scalar field, 0, that has a nontrivial 
dependence on the fifth dimension of the bulk spacetime. The gravitational action is 





(4) 
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where the gravitational constant k 2 is a free parameter ~ 0(1/N%) and V(4>) defines 
the nontrivial profile displayed by <fi in the bulk. Conformal invariance in the UV can 
be obtained when <f> — > at the boundary (u — > 0) and 

lim V(<f>) = ~ + ^A(A - 4)0 2 + 0( 4 ), (5) 

where the mass squared of the scalar m 2 i? 2 = A(A — 4) > —4 according to the 
Breitenlohner-Freedman (BF) bound [12J. Here A is the UV scaling dimension of the 
operator that couples with <\> at the boundary Eq. ^ implies that the spacetime will 
be asymptotically AdS 5 with radius R. 

The most general ansatz for the metric is [H] 

ds 2 = e 2A(«) ( _ /(u)df 2 + rff2) + e 2B{u)^_ ^ = ^ u y (g) 

J\ u ) 

(e B has dimensions of length). The boundary of the bulk spacetime is located at u — 0. 
Finite temperature effects are included by considering solutions of the classical equations 
of motion from Eq. Q with a horizon at u = Uh where f(uh) = 0. The functions above 
are assumed to be finite at Uh- We use a gauge where = u [Hj. Using this ansatz for 
the metric the entropy density is given by the area of the black brane horizon 

3A(u h ) 



while the Hawking temperature is (we use a prime for d/du) 
T = e AK)-BK)| fK) | 

47T ^ ' 

The class of V(4>) potentials considered here is of the form [11] 



-12 cosli70 + 6i(/> 2 
V VP) = ^2 ' ( 9 ) 

where 7 is related to the speed of sound in the IR and hi determines A near the UV 



fixed point. The choice 7 = yl/2 and b\ = 1.95 gives A = 3.3784 and leads to a 
black hole solution with the thermodynamic properties shown in Fig. 1, which matches 
the lattice data for c 2 in SU(3) YM above T c , as computed by Boyd et al. p2]. Note 
that even though the gravity calculations performed here require that N c — > 00 to avoid 
loop corrections, a comparison of our results to lattice data is still meaningful because, 
apart from a trivial rescaling, the thermodynamic properties of pure glue at iV c = 3 
and iV c > 3 are very similar, as shown in [2]. Our calculations are done by computing 
numerically the entropy density using Eq. ([7]). The other thermodynamic quantities 
are obtained from the entropy density via the thermodynamic identities and they also 
match the available data really well above T c . 

3. Polyakov Loops and the Gauge/String Duality 



When N c — > 00 and the radius of the background spacetime is much larger than the 
string length, an infinitely massive excitation in the fundamental representation of 
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Figure 1. cf(T) in the gravity model (solid red line) compared to (an interpolation 
of) the lattice data (dotted line) from [13] . 



SU (N c ) in the gauge theory at finite temperature is dual to a classical string in the 
bulk that goes from a probe D-brane at the boundary and to the black brane horizon 
[TH [T5J US] • The endpoints of the strings are non-dynamical fundamental probes and 
they are equivalent to the infinitely heavy quark limit in quenched lattice gauge theory. 

The correlator of Polyakov loops in this approach is given by (£(x)£*(0)) = Z string , 
where Z string is the full string generating functional that includes a sum over all the 
string worldsheets V whose boundaries describe a QQ pair separated by a distance r 
in the gauge theory [IB] . Within this approximation, the string dynamics (in Euclidean 
space) is given by the classical Nambu-Goto (NG) action 

S NG (V) =T S f dV^det hf (10) 
Jv 

where T s = ^ is the string tension, hf = G^d a X»d b X u (a, b = 1, 2), X" = X»(t, a) 
is the embedding of the string in the 5-dimensional spacetime, and G% v is the background 
bulk metric in the string frame. We assume here that the metric in these frames are 
related as follows: G? = eV^^G^ [TH US]. 

The disconnected contribution to the correlator in theories such as Q can be 
computed as follows. Using the gauge choice X M = (t, x(u), 0, 0, u) where the string 
endpoints at the boundary are at x = r/2 and x = —r/2, the regularized free energy of 
a single (static) heavy quark (up to a T independent constant) above T c is 

du [y/M$j- VMo(«)J - T s ut A (11) 

where M(u) = \G s oq\G SjUU is a positive- defined function of u that is regular at the 
horizon. Close to the boundary this function is Mq{u) = lim^o M(u) = u A - 4 / (4 — A) 2 
as long as V obeys Eq. 

Instead of using Eq. (11), it is easier [19] to work with its derivative dFq/dT = 
(dFq/duh) duh/dT. In fact, 

^ = T s JW^h) = T s e^ Uh e A ^+ B ^ (12) 
duh 
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and one obtains that 

^ = r . e ^^,_*£_ (13 ) 

Since the speed of sound squared in the plasma is c 2 = dlnT/dlias we see that 
duh/d\ns = c 2 s duh/ 'dhxT '. Thus, one can rewrite the equation above explicitly in terms 

d ^^T s e^ e Aiu h)+B{ u h) sdu, 
dT T cj ds 

It is possible to show [TU] that in gravity duals governed by Eq. Q the following equation 
holds 



dF Q e V3^( T ) 1 

~dT - Ws V{u h {T))cl{T)- (15) 

This shows that dFq/dT < (as expected from thermodynamics). Since V(uh) is well 
defined at finite Uh, one can see that dFq/dT only diverges when c s —> 0. Moreover, 
for confining gravity duals the magnitude of the jump in i at T c is determined by how 
quickly c s — > at T — > T c + . In the case of a first-order transition it can be shown 
explicitly that the thermal width of heavy quarkonia at strong coupling derived in [2"U] 
vanishes below T c . Continuous transitions where I smoothly interpolates between zero 
and its high-T limit, such as the one that occur in QCD with dynamical fermions |21j, 
can also be obtained depending on the choice for V (</>). Because of our choice for V, the 
medium is always in the deconfined phase and t can be evaluated at any T using (15) 
and there is no jump at T c . 

The renormalized Polyakov loop in quenched SU(3) QCD was computed some years 
ago on the lattice in [22] and more recently in [23]. We calculated I by integrating Eq. 
(15) and, because of the regularization, Fq is determined only up to a constant whose 
value is chosen via a match to the lattice data at T = 1.03 T c (the first data point above 
T c in Ref. [23] for N t = 4) for a given choice of a'. Our results for the renormalized 
Polyakov loop are shown in Fig. 2 (with a' /R 2 = 3.5) and one can see they agree very 
well with the lattice calculations from [23] between T/T c = 1—3. At higher temperatures 
the curve starts to considerably deviates from the data. This may be due to the fact 
that the gauge theory considered here is not asymptotically free. 



4. Conclusions 



We have shown that it is possible to describe the current lattice data for the 
thermodynamics and also the expectation value of the Polyakov loop of quenched SU(3) 
QCD using gravity duals where conformal invariance is softly broken by the presence 
of a nontrivial scalar field. The match gets worse as the temperature is increased, 
which may be a consequence of the fact that the theories considered here flow to a 
nontrivial fixed point (conformal invariance is restored when T — > oo but the system is 
still strongly interacting). It would be interesting to study examples where the gauge 



Polyakov Loops in Strongly-Coupled Plasmas with Gravity Duals 6 

1.4 
1.2 
1.0 
Lr 0.8 
0.6 
0.4 

2 4 6 8 10 12 14 

T/T c 

Figure 2. The renormalized Polyakov loop (solid black line) for the gravity dual in 
Eqs. Q and ^ and the lattice data for quenched QCD with N t = 4 (red) and N t — 8 
(blue) from Ref. [53]. 

theory is confining [17] and also investigate the properties of the Polyakov loop and its 
correlator in this case. 

Note that in order to match the lattice data of the loop a' / R 2 > 1, which indicates 
that higher-order derivative corrections [21] to the gravity dual are probably going to 
be important in this type of analysis. These corrections were used in [25] to study the 
dependence of I with the shear viscosity to entropy density ratio, r]/s, in superconformal 
strongly-coupled gauge theories. It would be interesting to extend this analysis to non- 
conformal gauge theories as well. In this case, rj and I will be temperature dependent 
and one could obtain how 77 varies with £ at strong t'Hooft coupling (this was computed 
at weak coupling in [26J). 
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